o 

<N 

I— i 

o 
o 



> 

o 
o 

p 

o 
o 



X 



Effect of spatial inhomogeneity on the mapping between strongly interacting fermions 

and weakly interacting spins 

Vivian V. Franca^' ^'Q and Klaus Capelle^'^ 

^ Physikalisches Institut, Albert- Ludwigs Universitdt, Hermann-Herder-Str.3, Freiburg, Germany 

^Instituto de Fisica de Sao Carlos, Universidade de Sao Paulo, Sao Carlos, 13560-970 Sao Paulo, Brazil 

^ Centro de Ciencias Naturais e Humanas, Universidade Federal do ABC, Santo Andre, 09210-170 Sao Paulo, Brazil 

(Dated: July 5, 2010) 

A combined analytical and numerical study is performed of the mapping between strongly inter- 
acting fermions and weakly interacting spins, in the framework of the Hubbard, t-J and Heisen- 
berg models. While for spatially homogeneous models in the thermodynamic limit the mapping is 
thoroughly understood, we here focus on aspects that become relevant in spatially inhomogeneous 
situations, such as the effect of boundaries, impurities, super lattices and interfaces. We consider 
parameter regimes that are relevant for traditional applications of these models, such as electrons 
in cuprates and manganites, and for more recent applications to atoms in optical lattices. The 
rate of the mapping as a function of the interaction strength is determined from the Bethe-Ansatz 
for infinite systems and from numerical diagonalization for finite systems. We show analytically 
that if translational symmetry is broken through the presence of impurities, the mapping persists 
and is, in a certain sense, as local as possible, provided the spin-spin interaction between two sites 
of the Heisenberg model is calculated from the harmonic mean of the onsite Coulomb interaction 
on adjacent sites of the Hubbard model. Numerical calculations corroborate these findings also in 
interfaces and super lattices, where analytical calculations are more complicated. 

PACS numbers: 75.10.Jm, 71.10.Fd 



I. INTRODUCTION 

Strongly interacting fermions are part of some of to- 
days most studied physical systems. In cuprate and 
manganite systems, for example, strongly correlated elec- 
trons are held responsible for high-temperature supercon- 
ductivity and colossal magnetoresistance, respectively^"— 
Systems of strongly interacting fermionic atoms can be 
realized in optical lattices, and are currently under in- 
tense investigation due to the possibility to use them 
as quantum simulators for understanding phenomena of 
condensed matter physics. ^^ 

Since the seminal works of Wigner on the low- 
density electron crystal^ and of Mott on the metal- 
insulator transition^ it is known that strong repulsive 
particle-particle interactions suppress itineracy and favor 
localization^ In the localized state, the repulsive inter- 
action is minimized, and charge degrees of freedom are 
frozen out. The dominating interactions in this state are 
magnetic. 

Mathematically, strongly interacting fermions are fre- 
quently described by the Hubbard model, which in one 
dimension is defined by the Hamiltonian 



The leading term of this expansion is the t-J model. 
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where Si is the spin one-half vector operator at each site. 
This model is frequently taken to be the starting point 
in investigations of doped cuprates. 

For a half-filled system, in which the number of 
fermions N equals the number of lattice sites L, the aver- 
age density n = N/L is unity. Since there are no empty 
sites, hopping is suppressed, and the t-J model reduces 
to the antiferromagnetic Heisenberg model 
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where cL and Cin 



are fermionic creation and destruction 
operators, hia = c\^Cia- is the particle-density operator, 
U the onsite interaction and t the hopping parameter. 

For sufficiently strong interactions, the Hubbard model 
can be expanded in powers of t/U . (Below we quantify 
what interactions can be considered 'sufficiently strong'.) 



where J = At^ /U and charge fluctuations are completely 
frozen out. The original system of strongly interacting 
itinerant fermions {U /t ^ 1) has thus been mapped on 
a system of localized spins with weak antiferromagnetic 
interactions (0 < J/t <C 1). 

The mathematics and the physics of this mapping 
are very well understood and discussed in the textbook 
literature. ^^'^^ The mapping of strongly interacting itin- 
erant fermions on weakly interacting localized spins is a 
standard concept of condensed-matter physics, routinely 
used in the interpretation of experiments on strongly cor- 
related solids. Recently, however, three important classes 
of systems have been discovered or created that call for 



a reconsideration and more detailed investigation of this 
mapping. 

It is well known that many strongly correlated sys- 
tems are characterized by nanoscale spatial inhomogene- 
ity. Such inhomogeneity can take the form of irregu- 
lar spatial variations of system properties, such as ob- 
served by scanning-tunneling microscope techniques in 
many cuprates and similar materials, ^^^^^^ or the form 
of regular spatial variations such as in naturally occurring 
or man-made superlattice structures J^"— In the presence 
of either type of inhomogeneity, the parameters charac- 
terizing the model Hamiltonian become site dependent. 
In the simplest case, with which we are mostly concerned 
here, the above homogeneous Hubbard model is replaced 
by an inhomogeneous model of the form 
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in which the onsite interaction Ui varies from site to site. 

A second class of systems we are concerned with here 
are nanoscale devices. In the modeling of such devices 
inhomogeneities in the model parameters occur simply 
because on the nanoscale the effect of the surface can 
no longer be neglected, and also because a typical de- 
vice combines more than one material, with the resulting 
interface automatically implying the existence of spatial 
variations in the system parameters. 

Finally, in still another line of research, ultracold atom 
gases have been trapped optically and arranged in optical 
latticesj^i^ In optical traps the system parameters can be 
controlled and varied in ways not possible in solid-state 
situations. In particular, the onsite interaction Ui can 
attain values U/t ^ 100 or larger. Such values are way 
beyond what is considered 'strongly correlated' in solid- 
state physics. 

Motivated by all these systems, we present, in the 
present paper, a combined analytical and numerical 
study of the mapping from the Hubbard model to the 
Heisenberg model in the presence of spatial inhomogene- 
ity. In Sec. ini we investigate the rate of the mapping, 
as measured by the difference in ground-state energies of 
the Hamiltonians. We allow the interaction strength to 
go beyond its typical solid-state values and to enter the 
ultrastrong regime attainable for cold atoms. 

In Sec, mil we turn to our main subject, the inhomoge- 
neous Hubbard model of Eq. (|4]). In Sec. lIII Al we investi- 
gate analytically the case of a single impurity, described 
as one site with a value of U differing from all others, 
and show that the Hubbard-to-Heisenberg mapping is 
preserved essentially in its homogeneous form, provided 
the effective J is calculated from the harmonic mean of 
the values of U on the sites connected by J. We illus- 
trate this finding numerically, by contrasting, for an im- 
purity system, results obtained from the harmonic mean 
with results obtained from the arithmetic, geometric and 
quadratic means. In Sec. lIII Bl we show that the harmonic 
mean allows to extend the Hubbard-to-Heisenberg map- 
ping to systems with more complicated types of spatial 



FIG. 1: Per-site GS energy as a function of interaction in the 
thermodynamic limit. Comparison between the Heisenberg 
with J = At^ /U and the half- filled Hubbard chain. The main 
figure includes values of U that can be realized in systems of 
trapped cold atoms, while the inset displays results for values 
of U typical of weakly and strongly correlated solids. 



inhomogeneity, such as super lattices, disordered systems 
and interfaces between different materials. Section IIVI 
contains our conclusions. 



II. RATE OF THE MAPPING FOR 
TRANSLATIONALLY INVARIANT SYSTEMS 

In a first step, to provide the background for the later 
investigations, we consider spatially homogeneous infi- 
nite Hubbard and Heisenberg chains, and investigate the 
rate of the approach of the ground-state energies of both 
models as a function of U . This allows us to quantify the 
rate at which charge fluctuations are frozen out. 

The per-site ground-state (GS) energy of the Heisen- 
berg chain in the thermodynamic limit is 
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The per-site GS energy of infinite half- filled (n ■ 
bard chain at U/t ^ oo is 
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Both expressions become identical for J = At'^/U. In 
order to quantitatively investigate the mapping at finite 
U/t^ we calculate e^^^^(n = 1, t/) by numerically solving 
the Bethe-Ansatz integral equation^*^ as a function of 
U and compare the result to the energy of the Heisenberg 
model at the corresponding value of J, i.e. e^^**(J = 
At^/U). The result is displayed^! in Fig. [H 

In Table [J we show the relative percentage deviation 
between the GS energies. 
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TABLE I: Relative percentage deviation of the per-site GS 
energies for values of U typical of weakly correlated solids 
(U — 1), strongly correlated solids {U — Q) and values that 
are attainable for atoms in optical traps. 
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for various representative values of U. This compar- 
ison between both models becomes trivial, once the 
Bethe-Ansatz solution is available, but already leads to 
a first somewhat unexpected conclusion. Frequently, 
the Heisenberg model is taken to be the starting point 
for a description of undoped antiferromagnetic insulat- 
ing parent compounds of high-temperature superconduc- 
tors. The effect of doping is accounted for by going from 
the Heisenberg to the t-J model, arguing that the latter 
should be a reasonable approximation to the Hubbard 
model for the involved large values of U. 

What the comparison in Fig. [1] and Table |I] shows is 
that for values of U that are representative of cuprate 
materials the t-J or Heisenberg models provide at best a 
semiquantitative approximation to the Hubbard model. 
At U = 6 the difference between both ground-state ener- 
gies is approximately 10%. Charge fluctuations are thus 
not yet frozen out for such [/, even at half filling. The 
rather large deviation observed shows that the mapping 
of strongly interacting fermions onto weakly interacting 
spins is not quantitatively reliable for, e.g., cuprate sys- 
tems at realistic values of U. There is no doubt, of 
course, that the t-J model captures the correct physics 
of the large-U Hubbard model - the above questioning 
only refers to the accuracy to which one needs to ob- 
tain a solution of the former, given that in the parameter 
regime typical of strongly-correlated solids it is itself only 
a moderately accurate representation of the latter. 

We note that this analysis is based on GS energies. 
An alternative comparison between both Hamiltonians 
would proceed in terms of the overlap of their wave func- 
tions, instead of the difference of their energies. Our 
main interest in this initial investigation is to investigate 
for which values of U the mapping breaks down, and for 
this it is enough to find one quantity that is not prop- 
erly reproduced. Thus, if we use our analysis to indicate 
when the mapping does not hold, we are on the safe side 
by using energies. Still another, and indeed more funda- 
mental, mode of analysis proceeds by directly comparing 
the Hamiltonians. We use this procedure in Sec. IHI A\ 
where our interest is not only in when the mapping breaks 
down, but also in how it can be restored. 



III. MAPPING IN THE PRESENCE OF 
SPATIAL INHOMOGENEITY 

We now turn to systems where the inhomogeneity oc- 
curs not only at the surface, due to finite size, but also 
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FIG. 2: Illustration of the inhomogeneous Hubbard model 
with a single impurity U' ^ U ^ and its putative relation to 
an inhomogeneous Heisenberg model with two bond defects 
J' ^J. 



in the bulk. A typical case is that of a localized impurity 
or defect, modeled by one site with onsite interaction dif- 
fering from that of all the others. We take this simple 
system as representative of Hubbard models with broken 
translational symmetry, and focus most of our analysis 
on it. The extension of our conclusions to interfaces and 
superlattices is discussed in Sec. IIIIBi 

Intuitively, one would expect that a localized pertur- 
bation of the homogeneous Hubbard model should only 
produce a similarly localized perturbation of the homoge- 
neous Heisenberg model. However, the relation between 
both models involves a projection on the subspace with 
no double occupation, together with an expansion in in- 
verse powers of f/^iii and it is not clear from the outset 
to what extent these operations preserve the above naive 
expectation of locality. 

In fact, there is one sense in which the mapping, if 
it continues to exist, cannot be local: U is defined on 
one site, while J connects two adjacent sites. While this 
difference is almost irrelevant in the homogeneous case 
where all sites are equivalent and translational symmetry 
rules, it becomes important in the inhomogeneous case, 
where any change in the onsite U on the Hubbard model 
must affect the corresponding intersite J for at least two 
sites of the Heisenberg model. This is illustrated in Fig.[2l 

The questions to address are thus (i) whether the map- 
ping still exists in the absence of translational symmetry, 
(ii) how to calculate the Heisenberg J from the Hubbard 
U in inhomogeneous systems, and (iii) if the mapping 
is as local as possible, i.e., involves only sites adjacent 
to the impurity site, or requires a higher degree of non- 
locality. In next subsections we address these questions 
analytical and numerically. 



A. A single impurity 

We start with the one-dimensional Hubbard model in 
the presence of a single impurity at site k with onsite 
interaction U' differing from the background value U on 
all other sites i ^ k^ 
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The standard proof of the mapping from the Hubbard to 
the t-J modeliSiii can be repeated for this Hamiltonian 
and leads to the inhomogeneous t-J model with Hamil- 
tonian 
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where we assumed that both U and U' are much larger 
than t. We rewrite this Hamiltonian by extracting one 
term from the sum over i ^ k io obtain 
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Now we take the average density n = N/L = 1. A priori 
this average can be obtained from many different dis- 
tributions Ui. However, since we are already in the limit 
[/, U' ^ t, the total interaction energy on the background 
and the impurity sites is minimized by the particular dis- 
tribution Hi = 1 for all i. Deviations from this are due to 
hopping processes that become increasingly suppressed 
as U and U^ grow. Thus the Hamiltonian ([TT]) at n = 1 
reduces to 
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where we chose / as a neighbor site of /c, k = I -\-l, such 
that Si • Si-^i = Sk • Sk-\-i' This, in turn, can be written 
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This derivation tells us that the Hubbard model with a 
single impurity can indeed be mapped onto a Heisenberg 
model with two bond defects, provided the impurity site 
and the background sites both have repulsive interactions 
that are much larger than t, and that the J connecting 
the impurity site with its neighbors to the left and to the 
right is calculated from the harmonic mean of the two 
onsite interactions. The mapping then still exists and 
is seen to be as local as possible, in the sense explained 
above. 

In order to investigate the rate of the mapping in inho- 
mogeneous systems, we now perform a numerical inves- 
tigation of both models. For illustration we also include 
in these calculations the quadratic, arithmetic and geo- 
metric means. 
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which has the form of a spatiahy inhomogeneous Heisen- 
berg model with background interaction J = At'^/U and 
two bond defects J' = At^ /U^ , where tj^ is the har- 
monic mean of U and [/', 
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Both the inhomogeneous Hubbard model with n = 1 and 
one U' ^ U and the inhomogeneous Heisenberg model 
with two J' ^ J are diagonalized numerically, and com- 
pared by means of the deviation between their GS ener- 
gies. We use the same method of analysis as in Sec. |TT1 
this time, however, applied to impurity systems. 

Our key result is contained in Fig.[3l which displays the 
relative percentage deviation between the inhomogeneous 
Hubbard and the inhomogeneous Heisenberg models with 
J' = J{U) for each of the four averages. The solid line 
represents the corresponding deviation obtained for the 
impurity- free Hubbard and Heisenberg models, where U 
and J = At^ /U are the same across the system. 

Figure [3]-(a) demonstrates clearly that if the back- 
ground U is so small that the mapping is not quanti- 
tatively reliable even in the homogeneous system (this 
includes the values found in cuprates), then it is also 
not reliable in inhomogeneous systems, regardless of the 
choice made for relating the bond defect to the impu- 
rity size. On the other hand. Fig. [3]-(b) makes a differ- 
ent statement: once the background U is large enough 
to permit the basic Hubbard-to-Heisenberg mapping to 
function, all four averages lead to values of J' = 4£^ /IJ 
that are equal (H) or different (A,G,Q) than obtained in 
the homogeneous system. 

This is unequivocal numerical evidence that the map- 
ping of strongly interacting fermions on weakly interact- 
ing spins survives in the presence of impurities and de- 
fects. In order to probe the locality of the mapping we 
have also performed numerical experiments with aver- 
ages over more than two neighboring sites, calculating 
J'{U) from, e.g. a weighted average of the interactions 
at the sites connected by U and their nearest-neighbor 
sites. No improvement (and frequently even worse re- 
sults) with regard to the simple two-site averages was 
obtained, indicating that the mapping is indeed local. 
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FIG. 3: Relative percentage deviation between the GS ener- 
gies of the single- impurity Hubbard chain and the two-defect 
Heisenberg chain, obtained from numerical diagonalization, 
adopting J' — J{U) with four different choices for the aver- 
age U: (a) values of U that are typical of solids and (b) val- 
ues of U that are attainable in systems of optically trapped 
atoms. System parameters: L = 8 sites, N = 8 fermions, 
open boundary conditions, impurity strength U^ = 3U/2. 



At first sight more surprising is that the alternative av- 
eraging procedures produce even smaller deviations for 
some values of [/, Fig. [3]-(b). However, for still larger 
values of U and U' (where the mapping should as a mat- 
ter of principle get better and better) all these alterna- 
tive averages produce deviations that continue to grow 
and yield D > 0, while the harmonic mean correctly ap- 
proaches the limit D = as /7, /7' -^ oo. This shows 
that only the harmonic mean has a chance to correctly 
describe the fermion-to-spins mapping in inhomogeneous 
systems, while the lower deviations of the other aver- 
ages for some parameter regimes only occur because the 
curves are monotonous so that there is always a range of 
values for which they are close to zero. 

We note that in Fig. [3] the ratio of U and U' was held 
fixed, such as to guarantee that for all values of [/, U^ 
was always substantially different from U {U' = 3U/2). 
In Fig.HJwe present the complementary analysis in which 
the background U is held fixed and the impurity inter- 
action is varied from U' < U io U' > U . For this com- 
parison we consistently adopted the harmonic mean. A 
first feature that jumps to the eye is that a single site 
with U' < U is enough to substantially deteriorate the 
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FIG. 4: Deviation of GS energies of the half-filled Hubbard 
model and the corresponding Heisenberg model as a function 
of impurity strength, obtained from the harmonic mean. At 
U' — U — Q the system is homogeneous (except for the pres- 
ence of the boundary) and the deviation is the same obtained 
in Sec. M 



mapping. By contrast, a single site with U' > U leads 
only to a slight reduction of the deviation between GS 
energies, much less than the deterioration observed for 
U' < U. 

This behavior arises from the hopping terms. Both at 
U' > U (more repulsive impurity site) and dX U' < U 
(less repulsive site) the on-site density at the impurity 
site slightly deviates from that at the background sites, 
as long as U and U' are both finite. In the latter case, 
however, hopping processes involving the impurity site 
increase as U' is reduced, and the Hubbard-to-Heisenberg 
mapping becomes correspondingly worse, while in the 
former case hopping continues to be strongly suppressed. 
The behavior displayed by Fig. |4]is thus consistent with 
what one would expect on the basis of the derivation 
leading from Eq.© to Eq.([l2l). 

Independently of, but in agreement with, our previ- 
ous analytical derivation we thus find that the harmonic 
mean solves the problem exactly for a single impurity and 
for sufficiently large interactions, while the other possible 
averages do not. However there is still one question con- 
cerning this result and it is addressed in next section: is 
the harmonic mean able to recover the mapping in more 
complex inhomogeneities? 



B. More complex inhomogeneities 

In view of our initial discussion of naturally occurring 
or man-made inhomogeneity in strongly correlated sys- 
tems, it becomes important to extend our analysis to 
more complex inhomogeneities than boundaries or single 
impurities. We here briefly describe our findings on three 
of these: interfaces, superlattices and disordered systems. 

Interfaces and superlattices can be described in the 
Hubbard and Heisenberg models as shown schematically 
in Fig. [5l While the description of superlattices by 



Hubbard Heisenberg 
interface /"V'W^V^'V'V^ 


U' U' U' U' U U U U ^ ^ 


J' J* J J* J' J* J 
superlattice /"V^^V^^^^'^V^Y^ 


U' U' U U U' U' U U ^ % 





FIG. 5: Illustration of simple interface and superlattice struc- 
tures in the Hubbard model, and their counterparts in the 
Heisenberg model. 



means of periodic spatial variations of U is the standard 
choice, ^^ which we here also adopt, it has been pointed 
out that a periodic modulation of local electric poten- 
tials can bring about a much larger change in the system 
properties ^ii^i2^ Here, however, we are interested in the 
Hubbard-to- Heisenberg transition, which is driven by the 
interaction and not by local potentials, and therefore we 
follow the usual prescription to ignore possible local elec- 
tric fields in the superlattice structure. 

We note that in the superlattice and the interface ge- 
ometry we now have three different spin-spin interactions, 
one, J, being calculated from /7, another, J', from U' 
and the last, J*, from U and U\ as indicated in Fig. [5l 
Figure [6] shows that for large U and U' the GS ener- 
gies of the Heisenberg model, when calculated from the 
harmonic mean, become identical to those of the corre- 
sponding Hubbard model, for all investigated geometries. 
In this sense, the mapping continues to work and to be 
as local as possible. Note that this would not be true for 
the arithmetic, geometric and quadratic means, whose 
deviations increase for large U and /7^ 

The fact that the Hubbard-Heisenberg deviation for 
the superlattice is larger than that for the interface can 
be understood on purely geometric grounds, as a conse- 
quence of the locality of the mapping and the nature of 
the harmonic mean: By comparing the distribution of J 
values. Fig. [5l we see that in going from the interface to 
the superlattice the number of interactions J and J' is re- 
duced by the same amount, while that of interactions J* 
increases. Since J' < J* < J, a reduction of the number 
of sites with J' worsens the mapping while a reduction of 
the number of sites with J improves it. To see what the 
net effect is we must take into account the interaction J* , 
which replaces J and J'. This interaction is calculated 
from /7, the harmonic mean of U and U' . The harmonic 
mean of any two positive numbers is less or equal their 
arithmetic mean, so that U is closer to U than to U^ and 
J* is closer to J than to J^ The substitution of an equal 
number of J and J' by J* thus has the effect of effec- 
tively increasing the number of 'bad' sites, and therefore 
to deteriorate the quality of the mapping. This is what 
the data show: the deviation for the superlattice is larger 
than that for the interfaces, if all other parameters are 
chosen the same. 
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FIG. 6: Relative percentage deviation between the GS ener- 
gies of the Hubbard and the Heisenberg model for a super- 
lattice structure, an interface, a single-impurity system and 
a system that is spatially homogeneous (except for the sur- 
face). In the superlattice and the interface system the num- 
ber of sites with U and U' is the same, respectively, the only 
difference being in their geometric distribution. System pa- 
rameters: L = S,U' = 3U/2, J' = At^/U' and J* = At^/tJ, 
where U is calculated from the harmonic mean of U and U' . 



This analysis shows that the harmonic-mean prescrip- 
tion continues to be usable for these more complex ge- 
ometries and that the effect of the geometrical distri- 
bution of U and U' sites across the system can be un- 
derstood and analyzed essentially on a site-by-site basis. 
This is a direct consequence of the locality of the map- 
ping. 

Disordered systems can be modeled simply by consid- 
ering a random distribution of impurities instead of just 
one. While a complete analysis of disordered systems re- 
quires statistical analysis of data resulting from a large 
number of realizations of the disorder, an analysis of a 
few representative cases is enough to conclude that the 
single-impurity results are not changed, in their essence, 
when the impurity concentration is increased. 

Specifically, we find that if all impurities have U' > U 
a higher concentration of impurities reduces the devia- 
tion between the GS energies. Keeping the concentration 
fixed and increasing U' /U also reduces the deviation, but 
to a much smaller degree. On the other hand, if all impu- 
rities have U' <U the agreement is naturally worsened. 
However, in the U' <U case the decisive factor is not so 
much the concentration of impurities but their strength, 
as measured by U' /U . This inversion in the effect of con- 
centration and strength of the impurities can be under- 
stood on the basis of Fig. IH which shows that for a single 
impurity with U' <U the deviation rapidly increases as 
U' becomes more different from [/, while for U' > U \i 
decreases only very slowly and almost saturates as the 
impurity sites effectively drop out of the system. 

This discussion shows that it is possible to control the 
degree of the fermion-spin mapping by means of the in- 



troduction of a suitable concentration of impurities of 
suitable strength. This possibility may be useful in the 
design of nanoscale devices based on strongly correlated 
systems, whose properties can be tailored from electron- 
like to spin-like by introducing suitable disorder. We note 
in passing that this is a strong-interaction effect, com- 
pletely different from the itinerant-to-localized transition 
resulting from disorder in Anderson localization. 



IV. CONCLUSIONS 

In the homogeneous situation, in which all sites are 
equivalent, the mapping is characterized by the behavior 
of the system as a function of the onsite interaction U. 
Not unexpectedly, the Heisenberg model is found to be a 
good approximation to the t-J and the Hubbard model 
at n = 1. Somewhat more unexpected is that the Heisen- 
berg model is rather a bad approximation to the Hubbard 
model at n = 1 even for values of U that are considered 
strongly correlated in solid-state applications. Only at U 
near 20 has the deviation dropped to about one percent. 
The standard mapping thus only becomes quantitatively 
reliable for values of U that are hard to reach in the solid 
state, but have already been demonstrated in cold- atom 
systems. 

In inhomogeneous situations, translational symmetry 



is broken. An analytical calculation for the simple case 
of a single impurity suggests that the mapping can be 
preserved in terms of the harmonic mean. Moreover, 
in terms of this mean the mapping is as local as pos- 
sible, i.e., the value of the Heisenberg J between two 
sites is only determined from the value of the Hubbard 
U at these two sites. Numerical calculations illustrate 
and corroborate this finding. This is more than a mathe- 
matical, or numerical, result: it means that the physics of 
the mapping, i.e., the gradual freezing out of the charge- 
degrees of freedom and the localization arising from the 
concomitant suppression of double occupation, is essen- 
tially the same regardless of the geometry and the pres- 
ence or absence of translational symmetry. 

The harmonic-mean prescription can be used easily 
and reliably for a wide variety of spatial inhomogeneities. 
Once the basic mapping is understood, the harmonic- 
mean prescription allows one to interpret and even to 
predict the behavior of much more complicated systems, 
without going through detailed analytical or computa- 
tionally expensive numerical calculations. 
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